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ABSTRACT 

Context. Soft X-ray high resolution spectroscopy of obscured AGNs shows a complex soft X-ray spectrum dominated 
by emission lines of He and H-like transitions of elements from carbon to neon, as well as L-shell transitions due to iron 
ions. 

Aims. In this paper we characterize the XMM-Newton RGS spectrum of the Seyfert 1 galaxy NGC 4051 observed during 
a low flux state and infer the physical properties of the emitting and absorbing gas in the soft X-ray regime. 
Methods. X-ray high-resolution spectroscopy offers a powerful diagnostic tool because the observed spectral features 
strongly depend on the physical properties of matter (ionization parameter U, electron density n e , hydrogen column 
density Nh), which in turn are tightly related to the location and size of the X-ray emitting clouds. We carried out 
a phenomenological study to identify the atomic transitions detected in the spectra. This study suggests that the 
spectrum is dominated by emission from a photoionized plasma. Then we used the photoionization code Cloudy to 
produce synthetic models for the emission line component and the warm absorber observed during phases of high 
intrinsic luminosity. 

Results. The low state spectrum cannot be described by a single photoionization component. A multi- ionization phase 
gas with an ionization parameter in the range of log U ~ 0.63 — 1.90 and a column density log Nh = 22.10 — 22.72 cm -2 
is required, while the electron density n e remains unconstrained. A warm absorber medium is required by the fit with 
the parameters log U ~ 0.85, log Nh = 23.40 and log n e < 5. The model is consistent with an X-ray emitting region at 

a distance > 5 x 10 -2 pc from the central engine. 

Key words, galaxies: Seyfert - galaxies: individual (NGC 4051) - Techniques: spectroscopic 



1. Introduction 

It is commonly accepted that the center of active galax- 
ies (Active Galactic Nuclei -AGNs) hosts a massive black 
hole (with a mass in the range 10 6 — 10 9 M ) accreting 
the surrounding material via the formation of a disk. How 
the energy released from the central engine interacts with 
the local environment and contributes to the history of the 
host galaxy is one of the crucial question of present astro- 
physical research. In this respect, while the mechanisms of 
energy output in the form of radiation and relativistic jets 
are quite well understood, it also seems that the outflowing 
winds have an important role in the overall energy budget. 
Although the origin of these winds is still controversial, at 
our present level of understanding the narrow-line regions , 
the inner part of an obscuring to rus (jBlustin et alj 



egions, 
. 2005) 



and the black hole accretion disk ((ElvisJ, 2000) are all pos- 
sible locations. 

X-ray obscured AGNs (with an intrinsic column density 
Nh > 10 22 cm -2 ) are not completely dark in the soft X-ray 

band. High resolution XMM-Newton and Chandra obser- 
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vations revealed a complex spectrum dominated by emis- 
sion lines from He-and H-like transitions of elements from 
carbon to neo n as well as by L-shell transitions of Fe XVII 
to Fexxi ions (ISako et al 112000 al. iKinkhabwala et af1l2002l 
iSambruna et alT 20011 Armentrout et al.ll2007D . This gas, 
which shows the signatu r e of a photoionization pro cess 
(jKinkhabwala et al J 120021 IGuainazzi fc Bianchil I2007D . is 
sometimes referred to as a warm mirror. 

In unobscured AGNs a modification of the output en- 
ergy spectrum may also occur as a consequence of ab- 
sorption by a warm ionized gas along the line of sight. 
The properties of these so called warm absorbers can be 
summarized as follows: i) average ionization parameter 
in the range log£ = — 3, ii) total column density in 
the range log Nh = 21 — 22 cm " 2 , iii) outflow velo cities 
of h undred of km s~* (see e.g. iBlustin et al.l 120051 but 
also ISteenbrugge et al.l I20091 ). Evidence of a multi-phase 
warm absorber gas was also recently reported for Mrk 841 
()Longinotti et al.ll20Toh . 

In general, detecting warm mirror signatures is easier in 
sources in low flux states, because the emission features are 
not outshone by the continuum radiation. This was the case 
for the Seyfert 1 galaxy Mrk 335, whose soft X-ray spectrum 
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resembled the spectra of obscured AGNs when the source 
was observed at low state (jLonginotti et aLll2008l) , but does 
not sh ow any evidence of a w arm absorber in the high flux 
state (|Longinotti et al.ll2007t ). 

The overall properties of the warm mirror (even if it is 
poorly constrained) and the warm absorber (as described 
above) are similar so that there is the possibility that they 
represent the same physical system. Conversely, the inter- 
play between the warm absorber and warm mirror regimes 
is best studied in sources that display both components. 

The source NGC 4051, a narrow- line Seyfert galaxy at 
the redshift of 0.00234, was at the center of many past in- 
vestigations in the X-ray band because it offers a unique 
laboratory where to test present theories and models about 
the physics of AGNs. T he X-ray em i ssion i s character- 
ized by rapid variations (|Lamer et all 2003; P onti et all 
2006 ) sometimes sho wing periods of low a ctivity (see 
Lawrence et all Il987l and lUttlev et all Il999tl . Its power 
spectral distribution (PSD) in high stat e resembles the 
beha vior of a galactic black hole system (|McHardv et al.l 
12004 . At high X-ray flux, the spectrum of the galaxy is 
characterized by a power law with photon index T ~ 1.8 — 2 
which becomes harder above 7 keV where a reflection com- 
ponent from cold matter has been observed. On long time- 
scales, the X-ray light curve of NGC 4051 shows low state 
flux periods of several months during which the spectrum 
in the energy range 2-10 keV becomes har der (r ~ 1) and 
show s a strong iron Ka line (as found by iGuainazzi et al.l 
1 19981 in Beppo-SAX data) . A soft X-ray excess is also evi- 
dent. 

As reported bv lQgle et al.l (|2004h . the high state X-ray 
spectrum of NGC 4051 in the soft band is a combination of 
continuum and emission line components. Curvature in the 
spectrum cannot be explained with simple models, i.e. a sin- 
gle power law or a black body, because an ionized absorber- 
emitt er has to be taken int o account as well. In this con- 
text, iKroneold et al.l (|2007h showed that the evolution in 
time of the properties of the warm absorber can constrain 
the physical parameters of the absorbing gas. In particu- 
lar they find that at least two different ionization compo- 
nents are required with matter densities of ~ 10 6 cm -3 and 
> 10 7 cm~ 3 , thus placing the warm absorber in the vicinity 

of the accretion-disk. Dynamical arguments permit us to 
infer that the warm absorber gas originates in a radiation- 
driven high-velocity o utflow in accretion disk instabilities 
(|Krongold et alJl2007h . 

On the other hand, as shown bv iPounds et al.l (|2004 . 
the low state flux spectrum of NGC 4051 is dominated by 
narrow emission lines and radiative recombination continua 
(RRC) from hydrogenic and He-like carbon, oxygen, neon 
and nitrogen. To be specific, a fit to the identified RRCs 
yields a mean temperature for the emitting gas of T ~ 4 x 
10 4 K, which favors a scenario invoking a photoionization 
process. In this case, the soft X-ray spectrum of NFG 4051 
in low state is similar to that o bserved for the prototyp e 
Seyfert 2 galaxy NGC 1068 (see lKinkhabwala et al J 120021 

Below w e do not repea t the a nalysis of the EPIC data 
but refer to IPounds et al.l (|2004 for more details on the 
main results obtained in the energy band 0.3-10 keV. We 
only say that a comparison between the EPIC PN data for 
the 2001 and 2002 observations shows that the high state 
observation flux level is a factor ~ 5 greater with respect to 
the low state. Furthermore, the spectrum shows a gradual 



flattening of the continuum slope from 3 keV up to 6.4 keV. 
It was also noted that when the fit to the 0.3-10 keV band 
continuum is extrapolated down in the soft X-ray (0.3-3 
keV) a strong excess appears in both the two observations, 
and as is clear from the RGS spectrum, it can be explained 
by a blending effect of fine structures (emission lines). 

Here we first conducted a phcnomcnological study of 
the emission lines identified in the spectrum of NGC 4051 
and compare our results with those known in literature. 
We further compared the RGS emission line spectrum with 
synthetic spectra generated w ith the photoionization code 
Cloudy 8 (jFerland et allll998l). For this purpose w e followed 
a similar approach as in lArmentrout et al.l (|2007l) (to which 
we refer for more details) on NGC 4151. 

The paper is structured as follows: in Sect. 2 we briefly 
describe the reduction of the XMM-Newton data set and 
describe our phenomenological analysis of the soft X-ray 
spectrum of NGC 4051. In Sects. 3 and 4 we give details on 
the Cloudy model developed and address some conclusions. 

2. A phenomenological study of the low state of 
NGC 4051: data reduction and line identification 

The source NGC 4051 {a = 12' i 03 m 09.6 s and 5 = 
44 d 31 m 53.0 s ) was observed by the XMM-Newton satel- 
lite on two occasions: on May 2001 for ~ 122 ks and on 
November 2002 for ~ 52 ks. While the former observa- 
tion coincided with a period where the central engine was 
brig ht (with luminosity of 7 x 10 41 erg s -1 in the 0.3 — 10 
keV. IPounds et alll2004 . the latter corresponded to a low 
state X-ray flux (corresponding to a luminosity of 1.5 x 10 41 
erg s^ 1 in the 0.3 - 10 keV, IPounds et all 12004 due to a 
low nuclear activity. This observation was conducted ~ 2 
days after the onset of the low state ([Pounds et al.l [20041 
Below we focus on the low state data analysis, because the 
warm absorber observed a t high state was a l ready well stud- 
ied with physical models (lOgle et aT]|2004 IKroneold et all 
l2007llSteenbrugge et al.l2009h . On the contrary, no attempt 
has yet been made to model the warm mirror in the low 
state with a self consistent physical model. 

The ODF files (OM, MOS, PN and RGS) were pro- 
cessed with the XMM-Science Analysis System (SAS ver- 
sion 8). Hence the raw data were reduced using SAS tasks 
with standard settings and the most update calibration files 
to produce the source and background spectra as well as the 
corresponding response matrices for the RGS cameras. 

We used XSPEC 12.5.1 (|Arnaud et all [20071) for our 
quantitative analysis and adopted the cosmological param- 
eters H = 70 km s _1 Mpc" 1 , Q, A = 0.73 and Q m = 0.27. 

To study the soft X-ray spectra of NGC 4051 in more 
detail, we then examined the first order spectra obtained 
by the XMM-Newton gratings. The spectral resolution of 
RGS in the first order spectrum is FWHM=72 mA and 
the calibration in wavelength is accurate up to 8 mA corre- 
spon ding to FWHM~ 620 km s" 1 and ~ 69 km s -1 at 
35 A ([XMM-Newton Users Handbook! I2009T ). Below we use 
the unbinned RGS 1 and RGS 2 spectra for the quantita- 
tive analysis. In Fig. [TJ we show the fluxed RGS spectrum 
in the wavelength range 5-38 A. 

As can be clearly seen, the 2002 RGS spectrum of 
NGC 4051 shows an unresolved continuum with a pre- 
dominance of emission lines (see for c omparison the RGS 
spectrum of NGC 1068 presented in iKinkhabwala et all 
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Fig. 1. Fluxed RGS spectrum of NGC 4051 (low state). The first order spectra of the two RGS cameras were combined 
and the resulting spectrum smoothed with a triangular kernel. The identified lines are labeled with the corresponding 
ion transition name and vertical dashed lines (the big dips at ~ 13 Aand ~ 21 Ain this plot are due to CCD gaps). 



2002). By contrast, the high state RGS spectra of the 
same source show a higher continuum flux level with a 
pronounced c urvature around ~ 1 5 Aand several absorp- 
tion features (jPounds et all . 120041) typical of a warm ab- 
sorber. Interestingly, the Nvi, Ovn and Neix forbidden 
lines are seen in bot h observations with the same flux level 
([Pounds et alJl2004h . 

The phenomenological spectral analysis follows the lo- 
cal fits method described in lGuainazzi fc Bianchil(|2007t ). In 
particular, the unbinned spectra are divided in intervals of 
~ 100 channels wide and Gaussian profiles are used to ac- 
count for all identified emission lines, with the line centroid 
energy as the only free parameter of the fit. 

Analogously, free-bound transitions (i.e. Radiative 
Recombination Continua) were also modeled as Gaussian 
profiles with free line width. Best-fit values of these widths 
are reported in Tabled] together with their errors. The local 
continuum was modeled as a power law with a fixed pho- 
ton index T = 1 and free normalization. For line triplets 
and for emission lines close to free-bound transitions, the 
relative distance between the central energies was frozen to 
the value predicted by atomic physics. 

We used the C-statisti c as the es timator of the good- 
ness of the performed fit (|Cashlll979ft . For any line under 
investigation, the emission feature was considered detected 
if, when repeating the fit without any Gaussian profile, we 
obtained a value of the C-statistic which differed from the 
previous one by at least 2.3, corresponding to the 68% con- 
fidence le vel (or, equiva l ently, 1 a for one interesting pa- 
rameter) (jArnaud et all 120071) . 

The results of the phenomenological fit to the emission 
lines are reported in Table [TJ We recall that the fluxes 
of the identified emission lines were estimated integrating 
over the Gaussian line profiles with which the emission fea- 



tures are modeled. This operation resulted in flux values 
somehow hi gher (up to a fact or of ~ 1.5) than those re- 
ported in iPounds et all (|2004l) even if the corresponding 
measured equivalent widths are fully consistent with the 
values quoted in the above mentioned paper. The origin of 
the discrepancy remains unknown. 



Here we give the best-fit parameters for the transitions 
identified in the soft X-ray spectrum of NGC 4051. For 
those lines that were not identified (i.e. AC < 2.3), but 
which are part still of triplet features, we give the upper 
limits to the corresponding flux value. We detected no ev- 
idence for either inflows or outflows, with an upper limit 
on the velo city of ~ 200 km s~ 1 , consistent with that es- 
timated bv lPounds et al.l (|2004). This was calculated from 
the full width half maximum of the distribution of residual 
velocity as derived from the difference between the expected 
and the measured laboratory wavelengths (see Fig. [2] and 
also Table [IJ. 



2.1. Radiative recombination continua (RRC) 



The electron temperature T e can be inferred by studying 
the profiles of the radiative recombination continua (RRC). 
In the RGS spectrum of NGC 4051 the RRCs detected with 
AC > 2.3 correspond to O VII, C v and C vi (see Table[2]for 
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Line ID 


^exp (A) 


A o6s (A) 


v (kms 1 ) 


Flux (xlO~ 14 cgs) 


AC 


Ne x Ly-a 
Neix (r) 
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-lo -,0+0.1)3 
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Neix (f) 
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21.600 


1, 7 +0-02 
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3-3±l;| 
7.0±i; 2 
< 0.80 


50 
21 
135 


Ovii (i) 
Ovii (f) 
N vn Ly-a 
Nvi (r) 
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22.101 
24.781 
28.787 


21.801°;^ 

22 10+ 03 
zz.iu_ 02 

24 7q+ U UJ 

28.78_ I ; 12 


1Q4+3400 
13y: -2380 
qo+400 
,3Z -250 
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13o -220 

-104+ 220 

1U ^-150 


i5-o±? : j 
Q+0.9 

z - u -0.8 

<0.92 


74 

238 

17 


Nvi (i) 
Nvi (f) 
Cvi Ly-a 


29.083 
29.534 
33.736 


29.07+001 
29 52+ - 02 

oo 7^+0.03 
' u _0.02 


_ q4 +220 
y4 -150 
_ 10 o+220 

9 qo+240 
^°' ;> _140 


1 q + 0-5 

L -°-0.6 

3-5±i;f 


6 

38 
29 



Table 1. Best-fit parameters for the transitions identified in the soft X-ray spectrum of NGC 4051. From left to right we 
give the name of the detected ion transition, the expected centroid wavelength (A) in the rest-frame as extracted from 
the CHIANTI database ((D ere et all I2001D . the observed wavelength (A), the corresponding shift in velocity (km s 1 ), 
the line flux (in units of 10 -14 erg s _1 cm -2 ) and the associated difference in C-statistics, i.e. AC. All lines that were 
not identified (AC < 2.3) but are part of triplet features are also reported. 



Line ID 




A o6s (A) 


v (kms 1 ) 


Flux (xlO -14 cgs) 


AC 


Ovii RRC 
Cv RRC 
Cvi RRC 


16.771 
31.622 
25.304 


16.78_^ 

O1.01_o 06 

25 10+ 010 


9 9n +8UU 
zzu -1200 

-wo±Z 

-1250±_gg 


3-8±i;« 
3-4±_;g 


38 
9 

24 



Table 2. The same as in Table Q] but for the identified RRCs. 



details) . Expr essing the temperature widths in eV to kT e 
(Liedhal 1999), we estimate them to be 

To vn = (5.0±t_) x 10 4 K , 

T c vi = (3.0±fj) x 10 4 K, (1) 
T c v = (l.2l.:|) x 10 4 K, 

respectively, so that the average gas temperature is T e = 

^3. 1+ 2 x 10 4 K , which agrees well with the result quoted 

by iPounds et alj (|2004h . It is to note that the low tem- 
perature found in this way is an indication that c ollisional 
ioniz ation and excitation processes are negligible (|Liedhall 
fl999h . 

In this phenomenological analysis, we used Gaussian 
profiles to fit the RRC features, which are in principle asym- 
metric. Still we verified that the use of a more appropriate 
model, as e.g. redge in XSPEC, gives consistent results. 

2.2. He-like triplet diagnostic 

We detected the most intense lines of He-like ions in the 
range 5-35 A. The transitions between the n = 2 shell 
and the n — 1 ground state shell as the resonance line 
(r : Is 2 1 S$-ls2p 1 -Pi), the two inter-combination lines 
(i : Is 21 So~ls2p 3 f_,i, often blended) and the forbidden 
line (f : Is 21 So-ls2s 3 Si) were detected. As demonstrated 



bv lPorauet fc Dubaul (|2000h the relative emission strength 
of the r, i and f lines are good indicators of the physical con- 
ditions of density and temperature of the gas. Using stan- 
dard notation we defined the ratios R = f /i, L = r/i and 
G={f + i)/r (jPorter fc Ferrand.ll2007lL Figure [3] shows the 
triplet of the O vn complex (forbidden, inter-combination 
and resonance lines) locally fitted by a power law and three 
Gaussian. In this case, following the phenomenological fit 
approach described in the previous Section, we only had a 
measurement for the fluxes of the / and i components (see 
Tabled]). 

With the flux measurements quoted in Table [TJ the 
previous relations give R = 2.lt° [l, L = 0.07 ± 0.06 
and G = 47^39 ■ Analogously, for the NVI triplet we get 
R = 2.7± 1 ; 2 ,, L<0.69 and G>5.4 (poorly constrained be- 
cause we only got an upper limit to the r line flux value), 
respectively. For Neix we had a lower limit only on the 
R(> 5.6) ratio, while the ratio of the O vni Ly-a to the 

Ovii forbidden intensity lines results in 0.4 7_n it- These 
line ra tios are consistent with the results bv lPounds et al.l 
(l2004h . 

2.3. Results of the phenomenological study: evidence of 
photoionized gas 

The results obtained from the phenomenological study al- 
low us to highlight some considerations on the physical con- 
ditions of the X-ray emi tting gas in NGC 4051. In deed, ac- 
cording to the study of iPorauet fc Dubaul ((2000) , a value 
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of the G ratio higher than 4 is a strong indication of a 
photoionized gas. An estimate of the gas electron density 
n e can be done when the other two line ratios L and R 
are taken into account. In the particular case of the Ovu 
triplet line ratios quoted above, the electron density is con- 
strained to be n e < 10 10 cm~ 3 for a pure photoionized gas 



(jPorouet fc D ubau 200(J ■ Note however that the line inten- 
sities obtained from the phenomenological study described 
above do not account for a warm absorber, which is not 
taken into account in the model. 

An additional constraint on the electron density value 
can be obtained noting that the XMM-Newton observation 
of NGC 4051 in its low state occurred ~ 20 days after 
the source entered this regime. Because the Ovu triplet 
line intensit y is consistent with that measured during high 
flux states ( Pounds et al.l l2004'). it is believed that the re- 
combination time of the O VII is larger than 20 days, thus 
implying (for a gas temperature of ~ 10 4 K) a more strin- 
gent constraint on the electron density of n e < 10 5 cm~ 3 

([Pounds et al.ll2004f ). 




Fig. 3. Zoom around the O vn triplet lines in the low flux 
state. The emission lines correspond to the f, i and r com- 
ponents, respectively. Only an upper limit to the resonance 
line r can be obtained (see text). The solid line represents 
the best-fit obtained with a power law + Gaussian lines 
model (details in text). 
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Fig. 2. Distribution of residual velocity for all lines in Table 
[1] with respect to the cosmological one. Note that the ob- 
served shifts are consistent with the cosmological ones, i.e. 
no outflow or inflow is observed. The solid line represents 
the Gaussian best-fit to the data. 



3. Fitting the spectra with the photoionization 
code Cloudy 

3.1. General properties of the model 

The results of the phenomenological analysis in Sect. 2 show 
that the bulk of the spectrum, measured by the RGS dur- 
ing the X-ray low state in NGC 4051, is dominated by 
phot oionization as alrea dy suggested by other authors (see 
e.g. lPounds et al.ll2004D . 

In this section w e use the photoionization code Cloudy 
(jFerland et al.l [19981) for modeling the overall spectrum of 
NGC 4051 in its low state assuming a plane parallel geom- 
etry with the central engine shining on the inner face of 



the cloud with a flux density depending on the ionization 
parameter U. 

The spectrum produced by a photoionized nebula crit- 
ically depends on the spectral energy distribution (SED) 
of th e ionizing cont i nuum . Below we adopt the AGN SED 
as in iKorista et al.l (|l997t) . In a typical AGN the observed 
continuum can be well represented by a SED characterized 
by several components: a big blue bump with temperature 
Tbb — 10 6 K (1 Ryd), a power law with a low energy expo- 
nential cut-off in the infrared region at kTiR = 0.01 Ryd; 
in the X-ray band (1.36 eV-100 keV) the SED is well ap- 
proximated by a power law with an exponential cut-off for 
energies lower than 1 Ryd; finally, for energies greater than 
100 keV an exponential fall as oc v~ 2 is usually assumed 
. We also included in the modeling the cosmic microwave 
background so that the incident continuum has a non-zero 
intensity for long wavelengths. Hence the AGN spectrum is 
described by the law 

F{i>) = iy auv e" p ^e~ ! ^ 1 + Av a "e~^^ , (2) 

where auv — —0.50, a x is the spectral photon in- 
dex and the constant A is obt ained requiring tha t 
F(2 keV)/F(2500 A) ~ 403.3"- (|Korista et al.l . Tl997h . 
where F(2 keV) and F(2500 A) are the flux densities at 
2 keV and 2500 A, respectively. 

To determine a x and u ox , we used the Epic data 
corresponding to the NGC 4051 high state observation. 
The resulting 0.2-10 keV energy band spectrum was fitted 
with a photoelectrically absorbed power law model within 
XSPEC, thus allowing us to measure a x — 0.96 ± 0.05 and 
F(2 keV) = (2.91 ± 0.01) x 10~ 29 erg s" 1 cm~ 2 Hz^ 1 , re- 
spectively. 

From the OM instrument we estimated the aperture 
photometry of the target in the UVM2 filter (centered at 
2310 A) obtaining a flux density of F(2310 A) = (1.50 ± 
0.01) x 10~ 14 erg s _1 cm~ 2 A -1 . From the flux densities at 
2 keV and 2500 A, the X-UV flux density ratio results in 
a ox ~ -1.14. 
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Once the AGN SED (erg s^ 1 cm~ 2 Hz -1 ) is known, it 
is straightforward to show that the number of hydrogen- 
ionizing photons the electron density n e and the di- 
mensionless ionization parameter U are related by U = 
Q/ Anr 2 n e c with r the distance between the central engine 
and the innermost illuminated layer of the clouds. Here we 
require that integrating over the SED (between 13.6 eV and 
13.6 ke V) we get the ioniz ing luminosity Li on ~ 4.1 x 10 42 
erg s -1 (|Qgle et all . I2004T ). In Fig. g] we compa re the SEP 
used i n this paper (solid line) with that given in lQgle et al.l 
(|2004l) . We recall that the dimensionless ionization param- 
eter U does not depend on the flux below 13.6 eV. The two 
spectral energy distributions give rise to comparable inte- 
grated fluxes in the 0.3-10 keV energy band (within a few 
percent). 

Note also that through the well known definitions of the 
ion izing luminosity Lj o n , of the number of ionizing photons 
Q (|Ferland et al.l fl998fl and the used SED F{ v), it is pos- 
sible to estimate a useful conversion relation between the 
dimensionless ionization parameter U and the ionizatio n 
parameter £ as given in iTarter. Tucker fc Salpeterl (|l969h . 
i.e. £ ~ 20U erg cm s _1 . 



-12 r 



-13 - 




log hi/ (keV) 

Fig. 4. SED used in the paper (solid line). Fo r compari- 
son w e also give (dashed line) the SED given in lOgle et all 
(2004). This was taken from the NED database and is likely 
to contain contamination from the hosting galaxy. 



3.2. Fit to the RGS spectrum 

Assuming the standard AGN continuum described above, 
we generated a grid of reflected spectra from a photoionized 
nebula, varying the ionization parameter log U, the elec- 
tron density logn e and the total column density log Nh- 
The free parameters spanned the ranges log U=[-l. 0,4.0]; 
log(n e /cm 3 ) = [2, 12] and log(A H /cm 2 )=[19,24] in steps of 
0.1 dex, respectively. 

Initially we extracted the line intensities from these 
simulations for all lines detected in the soft spectrum 
of NGC 4051. Hence, we tried to determine the best 
model that can describe the RGS line spectrum in the 
whole 5 — 35 Aband. Following the procedure described 

1 The default energy range used by the Cloudy code to eval- 
uate t he number of ionizin g photons Q is 1 Ryd - 7.354 x 10 6 
Ryd, (|Ferland et al.l . ll998D . 



in lLonginotti et al.l (|2008t ) , we calculated the value of the 
merit function for each grid model 

o 

where I is the intensity of each of the identified lines (with 
statistical error a a ) and I c is the intensity as predicted by 
Cloudy, both normalized to the value of the O vni Ly-a 
line. 

Minimizing the merit function quoted above gives a 
best-fit model (\ 2 ~ 13.8 with degrees of freedom v = 
10) corresponding to the parameter values \ogU = 0.4, 
log(n e /cm 3 ) = 4.4 and log(iV#/cm 2 ) = 21.8. A quantita- 
tive measure of the fit goodness for the used model is given 
b y the Chi-square Pr obability Function Q(x 2 , v) as defined 
in lPress et al.l (|2004t ). If the single phase component model 
is the true representation of the data, the probability to 
obtain the observed % 2 value is as high as Q ~ 12%. In this 
case, the model consisting of a single ionization state can 
be statistically rejected. 

We therefore investigated more complex models, includ- 
ing an additional warm mirror and one warm absorber cov- 
ering the combination of emitting components. For this ap- 
proach to be fruitful the constraints provided by the con- 
tinuum shape are crucial. Below we will fit the whole RGS 
spectrum globally. 

3.3. Global fit to the RGS spectrum 

We generated additive and multiplicative fits tables (with 
the same grid of parameters as before) to account for both 
the emission and absorption features observed in the RGS 
sp ectrum, and i mpor ted them within XSPEC as described 
in lPorter et alj (|2006f ). Our final model can be described by 
the formula phabs*mtab(n e , Nh, U)*{J2 atab(n e , Nh, U)}. 
Here, mtab and atab indicate the warm absorber component 
and the reflected component part of the spectrum depend- 
ing on the electron density, hydrogen column density and 
ionization parameter, respectively. In the model, the red- 
shift of each component is fixed to the cosmological value 
due to the lack of measurable velocity shifts from the phe- 
nomenological analysis (Fig. [2]) , while all the other parame- 
ters are free to vary. In the fit procedure we fixed the column 
density of neutral hydrogen to the average value observed 
in the Galaxy alo ng the line of sight to N GC 4051, i.e. 
1.32 x 10 20 cm" 2 (|Dickev fc Lockman|[T990T) . 

The fit does not formally depend on values of n e < 10 9 

cm~ 3 , which is expected because the ratios of the He- like 
triplets are insensitive to the electron density in th is region 
of the space parameter (|Porquet fc Dub au (2000)). Given 
the constraint on this parameter derived from the source 
time variability, we fixed its value to 10 5 cm -3 hereafter. 

We recursively increased the number of Cloudy additive 
components until this operation resulted in a statistically 
significant improvement of the fit quality. We found that 
two emission and one absorption components are required 
to fit the data. In particular, the final model corresponds 
to a value of the C-statistic of 6300 with 5178 d.o.f. and the 
model parameters are given in Table 4. Conversely, when 
the warm absorber component is not taken into account 
the fit visibly worsens and converges to a C-statistic value 
of 9452 with 5182 d.o.f. In this case, the line intensities 
corresponding to the He-like transitions are not correctly 
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estimated, with specifically the recombination line of the 
Ovn triplet well over-estimated (see Fig. [5]). 




21.8 22 
Wavelength (A) 



Fig. 5. Comparison of a model with two emission compo- 
nents (solid line) with the observed RGS spectrum (zoom 
around the Ovn triplet). 



We then extracted the line fluxes predicted by the best- 
fit model and compared them with the observations. In 
Figure |Hl we show with filled squares the intensities of 
all observed (see Table [lj and simulated (triangles) lines 
once normalized to the O vm Ly-a flux. The lower panel 
of the same figure shows the residuals between observation 
and theory. Note that the simulation underestimates the 
contribution of the Fexvn transitions, because the Cloudy 
database is in accurate for the corr esponding atomic param- 
eters (see e.g.[Bianchi et al.ll2010l ). The normalized intensi- 
ties of the observed lines as well as the Cloudy predictions 
are also reported in Table for clarity with the missing 
value of the Fe xvn transition. 

In Table 4 we give the relevant quantities estimated from 
the fit procedure (i.e. n e , Nh and U) together with their 
respective errors at the 90% confidence level for one inter- 
esting parameter. The 5-35 A low-state spectrum of NGC 
4051 is plotted in Fig. [7] with the best- fit Cloudy model su- 
perimposed on the observed RGS 1 (red) and RGS 2 (black) 
data and residuals in the lower part of each panel. 

Note that to estimate the covering factor of the source, 
we extracted the luminosity of the most prominent emis- 
sion line, Ovn (f), as predicted by Cloudy. For the two 
reflection components included in our model, we simulated 
the expected spectrum the SED described in Sect. 3.1 and 
fixing the electron density n e and hydrogen column den- 
sity Njj to the best-fit values given in Table 4. In addi- 
tion, the SED was normaliz ed to the i o nizing luminosity 
L lon ~ 4.1 x 10 42 erg s" 1 (|Qgle et all . I2004T ). Hence we 
fixed the distance from the source to the inner shell of the 
cloud to the lower values reported in Table 4 for each of the 
reflection components, i.e. d ~ 0.22 pc and d ~ 0.05 pc for 
the low and high component, respectively. After we defined 
the luminosity of the source, the e mission line lu minosities 
were predicted by Cloudy (see e.g.|Ferland 20081). 

With a redshift of 0.00234, the expected total inten- 
sity of the Ovn (f) line is ~ 10.6 x 10 3 erg s _1 cm~ 2 . 
Assuming a filling factor of unity, the ratio of the observed 
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Fig. 6. Intensities of the observed (filled squares) and simu- 
lated (triangles) lines corresponding to the best-fit and nor- 
malized to the O VIII Ly-a flux (see also Table [TJ. Residuals 
are in the right bottom panel. 



O VII (f ) intensity (see Table [T]) to the Cloudy expected 
value gives an estimate of the covering factor value, which 
turns out to be ~ 0.14. 



4. Discussion 

Most of the information on the physics and geometry of 
gas in AGNs is inferred by means of optical spectroscopy 
and imaging techniques with which it was shown that the 
AGN central high energy emission is the main source of 
ionizing photons with an occasional contribution from col- 
lisionally ionized plasma. In the last years X-rays obser- 
vations acquired an important role in AGN studies par- 
ticularly since Chandra showed the existence, at least for 
Seyfert 2 galaxies, of extended (a few kpc) X-ray emission 
(Bianchi et all 120061 ) similarly to what was observed in the 
optical band. High resolution spectroscopy in the soft X- 
ray band (0.2 — 2 keV) confirms the overall scenario, and 
photoionization seems to be the dominant ionization mech- 
anism which results in a spectrum characterized by recom- 
bination lines from He- and H-like transitions of C to Si 
elements and by Fe-L transitions. In this respect, X-ray 
high-resolution spectroscopy offers a powerful diagnostic 
tool because the observed spectral features strongly depend 
on the physical properties of matter (ionization parameter 
U, electron density n e , hydrogen column density Njj as well 
as size and location of the emitting clouds) . 

The Seyfert 1 object NGC 4051 shows a very rich emis- 
sion line X-ray spectrum when observed in low-flux state. 
According to the analysis we conducted on the XMM- 
Newton RGS data, the observed soft X-ray features orig- 
inate in a low-density photoionized gas. In order to con- 
strain the physical properties of the photoionized gas, 
we simulated synth etic spectra via the Cloudy software 
(jFerland et al.l[l998j) and compared them to the RGS data 
with standard minimization techniques. We found that to 
describe the overall soft X-ray spectrum, at least a three- 
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Line ID 


^exp (^-) 


Observed ratio 


Cloudy predicted ratio 


Ne x Ly-a 
Neix (r) 


12.134 
13.447 


2 +u ' 2 

u - z -0.2 

<0.1 


0.2 
0.1 


Neix (i) 


13.550 


<0.1 


0.2 


Neix (f) 
Fexvn 3s-2p 
viii Ly-a 
Ovii (r) 


13.700 
17.073 
18.969 
21.600 


u -'-0.2 

5 +0 ' 2 

U.0_ 2 

1 o +0 - 2 

<0.1 


0.6 

1.0 
0.3 


Ovii (i) 
Ovii (f) 
N VII Ly-a 
Nvi (r) 


21.790 
22.101 
24.781 
28.787 


i7) +o ' 2 

L - u -0.2 
9 9 +0.3 
^• z -0.3 

o-3±g : i 

<0.13 


0.6 
2.5 
0.2 
0.01 


Nvi (i) 
Nvi(f) 
C VI Ly-a 


29.083 
29.534 
33.736 


0.5t° ;» 


0.01 
0.1 
0.7 



Table 3. Intensities of the observed lines and of those simulated by Cloudy corresponding to the best-fit and normalized 
to the O VIII Ly-a flux (see text and Fig. [5] for more details). 





logU 


log(n e /cm J ) 


logiNn/cm- 1 ) 


d (pc) 


A/=i (pc) 


Low 
High 


63 +u us 

i.90t° o ;f o 


<5 
<5 


22 10 +U - 3U 

— 0.04 
99 ^9+0.25 

— 0.13 


>0.22 
>0.05 


>0.04 
>0.20 


Warm Abs. 


0.851^ 


<7 


23 36 +u ui 


>0.02 


>0.01 



Table 4. Best-fit parameters for the three components of the adopted Cloudy model used to fit the RGS data. Errors 
are given at the 90% confidence level for one interesting parameter. Columns are: the ionization parameter, electron 
density, Hydrogen column density, distance from the central ionizing source, and the average X-ray emitting source size 
evaluated as A ~ Nh/ fn e , where / is the filling factor. Here, for simplicity we have assumed / ~ 1. 



phase gas is required (two emission components and one 
warm absorbing component). Referring to the emission 
components respectively as low and high ionization com- 
ponents, our fit procedure gave us their physical proper- 
ties. For the low component we have \ogU ~ 0.63, and 
log(iVff/cm 2 ) ~ 22.10 and for the high component we 
have logC/ ~ 1.90, and log(iV ff /cm 2 ) ~ 22.20. Using 
Cloudy we get for the electron density n e an upper limit of 
log(rt e /cm 3 ) ~ 9, which reduces to log(rt e /cm 3 ) ~ 5 when 
the recombination time scale of O vn is taken into account. 
Even if the warm absorber gas seems to be required by our 
fit procedure, its parameters are poorly constrained. Thus 
it is characterized by log[/ ~ 0.85, log(7Vy /cm 2 ) ~ 23.36, 
and log(n e /cm 3 ) < 7. 

This technique was successfully applied before to the 
Seyfert 1 Mrk 335 and the Sey f ert 2/starburst galax y 
NGC 1365 (jLonginotti et all 120081: [Guainazzi et al.ll2009f) . 
The main difference is that NGC 4051 is characterized by 
a strong warm absorber component in the high flux state 
that is still affecting the spectrum even when the nuclear 
flux is attenuated. Indeed, we found out in our analysis of 
the low flux state data that the effect of the line of sight 
medium is not negligible, particularly not in the modeling 
of the resonance line of the O vn triplet (see Fig. [5]) which 
is close to several absorption features. For example, the 
resonance line cou ld be weakened by th e same line in ab- 
sorptiorQ (see e.g. Kroneol d et al.l [2007'l. Nonetheless, the 
physical parameters of the warm absorber cannot be well- 

2 As noted by iSako et all (|2000 bh and iKinkhabwala et alj 
(2002), the resonance line of the O vn triplet could be also en- 



constrained by the analysis of the low flux state data (see 
Table 4). 

The average distance r of each of the photoionized 
plasma-emitting components from the nuclear source can 
be estimated by the definition of the ionization parame- 
ter U after normalizing to the ionization luminosity Li on . 
However, our results are insensitive to values of the electron 
density n e lower than 10 5 cm -3 . In this limit, we can only 
determine a lower limit of the A-ray-emitting gas location 
(Table 4). 

The analysis carried out in this paper allowed us to 
identify two ionization states for the line emitting gas and 
one warm absorber medium. It is interesting to note that 

- The X-ray emitting region can be placed at a distance 
of r > 0.05 pc. 

Indeed. \6e\eeTaT\ rf200"i found that the NGC 4051 X- 
ray narrow-line regions can be placed at a distance of the 
same order of magnitude. This was also con firmed by the 
Chan dra ACIS-S images of the same galaxy ()Uttlev et al.l 
1999), which showed a size of the diffuse emission smaller 
than that of the optical narrow- line regions (30 -220 pc, 
iChristopoulou etaill997l) , thus implying a clear separation 
between the X-ray and optical emissions. 

This is also naturally exp ected as a consequen ce of 
projection effects: as shown by ISchmitt et al . (2003), who 

hanced by photoexcitation. Note however that this would also 
result in a boost of all the higher order resonance transitions 
of the H-like and He-like ions (Ly-/3, Ly-5, He-/3 and He-5), but 
this enhancement is not currently observed. 
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studied a sample of 60 Seyfert galaxies with the Hubble 
Space Telescope, the Seyfert 1 narrow-line regions objects 
are more circular and compact than those in the Seyfert 
2 galaxies, with the Seyfert 2 subsample characterized by 
more elongated shapes. This agrees well with the unified 
picture according to which the conical narrow-line region 
of a Seyfert 1 galaxy is observed close to the axis of sym- 
metry, while that of a Seyfert 2 galaxy is observed from an 
orthogonal line of sight. 

Furthermore, the scale-length found in this paper is 
consistent with t he inner radius of th e torus in NGC 4051 
as determined bv lBlustin et al.l (|2005l ). i.e. r ~ 0.15 pc. 

- The NGC 4051 low state warm absorber is poorly con- 
strained but its existence is nevertheless required by the fit. 
In particular, we found a lower limit of the warm absorber 
distance ~ 0.02 pc, i.e. at least a factor 10 larger than that 
measured in the high state flux (jKrongold et al.ll2007f) . 

Indeed, by using the long XMM-Newton exposure of 
NGC 4051 in its high flux state and studying the time 
evolution of the ioniz ation states of the X-ray absorbers, 
iKrongold et "all (l2007h were able to put severe constraints 
on the physical and geometrical properties of the warm ab- 
sorber medium. 

They specifically found that the warm absorber consists 
of two different ionization components which are located 
within 3.5 lt-days (or 0.0029 pc) from the central massive 
black hole. This result allowed the authors to exclude an 
origin in the dusty obscuring torus because the expected 
dust sublimation radius^] is at least one order of magnitude 
larger. Hence the authors suggested a model in which the 
black hole accretion disk is at the origin of a X-ray absorber 
wind, which forms a conical structure moving upward. 

If this is the correct picture, when the continuum source 
is switched off, the compact warm absorber might not be 
observed anymore during the low state flux of NGC 4051. 
Our analysis showed instead the existence of a more ex- 
terior X-ray absorber, which absorbs the soft X-ray pho- 
tons emitted from sources (as for exampl e the inner sur- 
face of the conical structure proposed by IKrongold et al.l 
|2007|) located (in projection) at scales larger than the torus 
and/or the narrow-line regions. Remarkably, this could in- 
dicate the existence of a diffuse warm material filling the 
wind-generated cone. 

Figure [5] gives a qualitative representation of the model. 
During the high state flux (left panel) a two ionization com- 
ponent warm absorber (here labeled as II) lying within 
a few 1-days (~ 0.003 pc) from the a c creting black hole 
was identified by IKrongold et al.l (|2007t ). lOgle et al.l (|2004T) 
found that the NGC 4051 X-ray narrow- line regions can 
be placed at a distance of r > 0.02 pc, while the opti- 
cal narrow-line regions are on the scale of tenth of parsec 
(|Christopoulou et al.lll997l ). During the low state flux (right 
panel), the interior warm absorber might not be observed 
anymore since the central engine is switched off. A more 
exterior warm absorber (labeled as I) could now absorb the 
X-ray photons emitted from sources located on the scale 
larger than the torus and/or the narrow-line regions. 
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3 The torus inner edge has to be at a distance larger than 
the d ust sublimation r adius r«u&. In the particular case of NGC 
4051. IKrongold et all [|2007D found r sub ~ 0.01 pc. 
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Fig. 7. The 5-35 A RGS 1 (red data points) and RGS 2 (black data points) spectra are shown together with the best-fit 
Cloudy model described in the text and corresponding to the physical parameters reported in Table 4. 




Fig. 8. The cartoon shows qualitatively the location of the X-ray emitting and absorbing material (see text for details). 



